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Abstract: Prepare a two-species BEC in a perfectly phase-mixed state. By applying Rabi 
drives, one can tune the range of wavelengths of phase-separating excitations that are dynamically 
unstable. Yada yada yada, this determines the characteristic sizes of the eventual spin domains. 
The trapping potential is neglected because it makes life hard, but of course this is a terrible 
approximation, and the results are of only inspirational value. Since this is thus a rather trivial 
calculation, the very modest increase in kudos that might accrue from journal publication does not 
seem to outweigh the grief of having to change the title. This will therefore be an ArXiv exclusive. 



/vq ' Consider a two-species condensate held in a homogeneous trap, and subjected to a Rabi drive. Homogeneous 

stationary states (there are always at least two solutions) satisfy the time-independent Gross-Pitaevskii equation 



A*ilk = Ae- m e ie ^ 2 + ( 5ll |</>i| 2 + g 12 1^ 2 | 2 ) Vi 

M2 t/- 2 = Ae^e-^Vi + (312 l^il 2 + 322 |V> 2 | 2 ) ^2 ■ (1) 

Here A is the driving amplitude, which initially we will take to be real and positive; Ol is the phase of the driving 
1 . ■ field; and tt is its frequency. (Setting A positive may obviously be done without loss of generality, but it will be 
£h \ notationally convenient later to redefine it so as to allow it to be negative.) 

Our goal is to determine the Bogoliubov spectrum for small perturbations around the homogeneous stationary 
states, and then to consider which Bogoliubov modes may be dynamically unstable. Such instabilities will indicate 
the earliest, linear stages of spontaneous phase separation into spin domains. By varying A, it will turn out to be 
rvq J possible to control the ranges of wavelengths that are unstable. If we assume that the nonlinear evolution of domain 
^ ■ formation does not significantly alter the 'decisions' made in the linear epoch, about where each species predominates, 
then control of the instability wavelengths means control of domain sizes Jy. (This is admittedly a yada theorem 
HH/; but it is probably not too far wrong.) 

From the viewpoint of BEC theory, it is obviously desirable to incorporate realistic features like an inhomogeneous 
trapping potential, and accurate values for the independent components of the interaction co-efficients gij . To get 
t— I I a feel for the basic concept, however, let us consider the tractable case of no potential, as in (|l|), and also stipulate 
<7n — 322 = fl- It will then be a nice notation to define 512 = <?A. We then have the two solutions 



C3 ; - t p 1= /p e *(0o+0z,/2) e -iGu+n/2)t 

^ 2 = ±^ e ^o-e L /2) e - l{ ^n/2)t 

I* = P(9 + 912) ±A= -Aj [1 + A/4 + a/2] , (2) 

g ■ 2Mi A 

where in the last step we introduce another notation that includes the healing length: g\ 2 = g(\ + A), A = ztagp, 
and 2gp = h/(2M£ i 2 ) for M the atomic mass. In these solutions the densities p of each of the two condensates are 
equal and constant, but their phases may differ by the drive phase, or by the drive phase plus tt. There are thus 
two stationary states to consider; and at least to a good first approximation, either of them may be prepared at will, 
with an appropriate initial Rabi pulse. (A caveat: it will ultimately be an important question how much dynamical 
instabilities may limit the accuracy of the initial state preparation.) So we have two Bogoliubov problems to solve, 
depending on the branch of the ± in (|2|). Since the positive parameter A already appears in our solutions, however, 
we have made things more compact by absorbing the ± into a, which can thus be positive or negative. So from 
now on we will discuss only one set of Bogoliubov modes, whose frequencies will depend on our arbitrary real control 
parameter a. 

If for j = 1 , 2 we write 

i>j - *i [l + e^ k ' x -^(r ik + s jk ) + e-« k — "^(r;,, - s* k )~ 
and then linearize the time-dependent Gross-Pitaevskii equation in r^ and Sjk, we obtain the Bogoliubov equations 
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This is trivially equivalent to 
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Hence we have a 'sonic branch' of excitations, in which the two species oscillate in phase, and in which the frequencies 
are always real (and in fact it can be shown, positive) as long as A is not too negative. For all condensates prepared 
so far, |A| <C 1, so these modes are not very interesting. 

On the other hand, we have an 'optical branch' in which the two species oscillate out of phase, and here is the 
physics we want to see. The frequency of an optical mode of wave number k is, for k = |k|, 



tu k = 2gp^/(k 2 e ~ ") (k 2 e 



A), 



(5) 



where the sign of any real uik does not concern us here. (A negative real frequency for a positive norm mode means 
only energetic instability, and we are interested in dynamical instabilities, for which the Bogoliubov norm is zero.) 
What matters to us is that u>k will be imaginary, indicating a dynamical instability, if 



a + A < fc 2 £ 2 < a , A < 
a < fc 2 £ 2 < a + A , A > 



(6) 



So with negative A, the two species intrinsically prefer to stay mixed; but by applying a positive a we can destabilize 
them toward phase separation. Positive A, on the other hand, means that the condensates would intrinsically prefer 
to separate, with domain sizes no smaller than £A -1 ' 2 (assuming that the minimum domain size is given by the 
maximum dynamically unstable fc). By applying negative a we can enlarge this minimum domain size; if we make 
it larger than the sample size, then we effectively frustrate the phase separation. 

For either sign of A, we can turn on phase separation by applying positive a. We can make the minimum domain 
size smaller; and for large enough a, we can even introduce a maximum domain size. Hence we can truly be masters 
of our domains. 



[1] The title refers to an ingenious Seinfeld exercise in evading censorship. 

[2] 'Yada yada yada' is a mysterious Seinfeld catchphrase which glosses over highly nontrivial details, relying on sheer speed 
to solicit the audience's belief in the conclusion. Hninim. 

[3] The assumption that the initial linear epoch controls the final domains is not obviously unreasonable; but it is not obviously 
valid, either. Even within Gross-Pitaevskii, the assumption needs to be checked by following the full, nonlinear evolution, 
with dissipation. Moreover, in the linear epoch there will typically be several competing unstable modes, and excitations 
of all of them will grow in quantum superposition, until mean field and Bogoliubov theory break down (after a time which 
grows only logarithmically in the particle number N). This then raises a deep and open problem, hitherto raised mainly 
in cosmology: in a nonlinear quantum field theory, how do quantum fluctuations at early times 'mature' into classical 
inhomogeneities at late times? So the experiment would be very interesting from a fundamental point of view. 



